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Analysis of the MOST light curve of the heavily spotted 
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ABSTRACT 

Continuous photometric observations of the visible component of the single-line, K2IV 
spectroscopic binary II Peg carried out by the MOST satellite during 31 consecutive 
days in 2008 have been analyzed. On top of spot-induced brightness modulation, eleven 
flares were detected of three distinct types characterized by different values of rise, 
decay and duration times. The flares showed a preference for occurrence at rotation 
phases when the most spotted hemisphere is directed to the observer, confirming 
previous similar reports. An attempt to detect a grazing primary minimum caused by 
the secondary component transiting in front of the visible star gave a negative result. 
The brightness variability caused by spots has been interpreted within a cold spot 
model. An assumption of differential rotation of the primary component gave a better 
fit to the light curve than a solid-body rotation model. 

Key words: stars: individual: II Peg, RS CVn-type, flare, star spots, rotation. 



I INTRODUCTION 

Studies of II Peg (HD 2240 85, Lalande 46 867) began with 
the spectroscopic analysis of lSanfordl (|l92ll ). He found II Peg 
to be a late-type (K2), single-line spectroscopic binary (SB1) 
and determined its first orbital elem ents. Well defined , reg- 
ular light variations were noticed bv lChugainovl l|l976l ) who 
explained them by rotation of the primary component of 

II Peg with a cold spot on its surface. He also observed flares 
and concluded that this is a BY Dra-type binary system. 
However, subsequent pho tom etric and sp ectroscopic obser- 
vations bv lRucinskil (|l977T ) and lVogtl (|l979h led these authors 

* based on data from the MOST satellite, a Canadian Space 
Agency mission, jointly operated by Dynacon Inc., the Univer- 
sity of Toronto Institute of Aerospace Studies, and the University 
of British Columbia, with the assistance of the University of Vi- 
enna. 
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to conclude that the star is more akin to RS CVn-type sys- 
tems, although with an invisible less-massive component; the 
distinguishing features of the BY Dra and RS CVn type bi- 
naries (with, respectively, dwarf and sub-giant components) 
were being defined at that time. Changes of e quivalent width 
of the H q line with orbital phase analyzed bv lBopp fc Noah] 
|l980al ) confirmed the RS CVn classification. 

Since the 1980's, II Peg was one of the most frequently 
observed RS CVn-type stars. A model of multiple spots 
was used for analysi s of the availab l e light -curve data sets 
for the first time by iBopp fc Noahl (|l980bl) . The most de- 
tailed and complete s tudy of II Peg wa s presen te d in a se- 
ries of four papers bv lBerdvugina et all (|l998al lbl [l999ai rbh. 
In the current paper we utilize most of the parameters de- 
rived in the first paper of this series which was based on 
high-resolution spectra used to define a high-quality ra- 
dial velocity orbit. In brief, the essential physical param- 
eters of the primary (visible) star were found to be: T e // = 
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4,600 ± 100 K, log g = 3.2 ± 0.2, [Fe/H] = -0.4 ± 0.1, 
v sini = 22.6 ± 0.5 km/s, _Ri = 3.4 ± 0.2 Kq, spectral type 
K2IV, with ephemeris for conjunction (visible star behind), 
T conj = 2, 449, 582.9268(48) + 6.724333(10) x E, where E is 
an integer number of orbits. 

From the analysis of TiO bands and simultaneous photo- 
metric observations, the fictitious, entirely unspotted visual 
magnitude of the primary star was estimated at a relatively 
bright level V u = 6.9; we return to this matter later in the 
paper as it affects the results of our spot modelling. The or- 
bital inclination was estimated at 60° ± 10°, leading to the 
primary mass Mi = 0.8 ±0.1 Mq and implying that the 
secondary star is probably a main-sequence, late-type dwarf 
(M0-M3V) with mass M 2 « 0.4 ± 0.1 M©. The presence of 
a white dwarf in this binary s ystem was previously excluded 
bv lUdalski fc Rucinskil (|l982h on the basis of ultraviolet ob- 
serv ations made by th e IUE sp acecraft. 

iBerdvugina et all |l998bl ) presented multi-epoch im- 
ages of the primary component, obtained by means of the 
Doppler imaging technique. They found that the spot dis- 
tribution and spot parameters obtained from the spectral 
analysis are in good accordance with th ose derived solely 
from an alysis of photometric observations. IBerdvugina et al.l 
( 1999b|) discussed the "flip-flop" phenomenon, i.e. a shift of 
the maximum spot-activity to the opposite side of the stel- 
lar surface. The authors also concluded that - because the 
largest active area tends to be located on the hemisphere 
facing the secondary star - this component may play an 
important role in the magnetic phenomena in the system. 

The current paper presents analysis of continuous ob- 
servations of II Peg conducted using the MOST satellite 
during 31 days in September and October 2008 (Section [2]), 
a circumstance which permitted us to address the following 
issues: (1) Study of frequency and orbital-phase localiza- 
tion/orientation of flares in the system (Section [3]); (2) A 
search for grazing eclipses caused by the secondary (Sec- 
tion!?]); (3) Determination of the differential rotation of the 
visible star as its minute signatures are better defined for a 
long observing run (Section [SJ. 



2 OBSERVATIONS AND DATA REDUCTION 

The optical system of the MOST satellite consists of a 
Rumak-Maksutov f/6 15 cm reflecting telescope. The cus- 
tom broad-band filter covers the spectral range of 380 - 
700 nm with effective wavelength falling close to Johnson's 
l/band. T he pre-launch charact eristics of the mission are de- 
scribed by IWalker et al.l d2003T) and th e initial post-launch 
performance bv lMatthews et all (|2004l ). 

II Peg was observed from 15th September to 16th Oc- 
tober 2008, in HJD = 2, 454, 725 - 2, 454, 756, during 439 
satellite orbits over 30.877 days. The individual exposures 
were 30 sec long. Only low stray-light orbital segments were 
used, lasting typically 25 min of the full 103 min satellite 
orbit. In spite of the high background, telemetry and South 
Atlantic Anomaly breaks, the almost continuous light curve 
is very well defined (Figure [1]). 

Because II Peg is usually close to or slightly fainter 
than 7th magnitude, i t was observed in t he direct-imaging 
mode of the satellite (| Walker et al.l 120031 1 . The CCD cam- 
era does not have a mechanical shutter which limits pos- 
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Figure 1. The light curve of II Peg in magnitudes. The horizontal 
scale is in heliocentric Julian Days (lower edge) and in orbital 
phase units (upper edge) with zero phase for conjunction with 
the visible star behind. The phases were calculated by means of 
the ephemeris determined bvfBcrdvugina et al . (1998a) as given 
in Section [l] 

sibilities of obtaining calibration frames, as it is com- 
monl y practised during ground-based observations. How- 
ever, iRowe et al.l (|2006al rbh proposed an excellent calibra- 
tion procedure: Because the background level caused by the 
Earth stray light usually changes very significantly during 
the orbital motion of the satellite, it is possible to deter- 
mine both the dark-level and the flat-field information for 
pixels within small images (rasters) around stars on a per- 
pixel basis. We removed first the background gradient vis- 
ible in most frames and caused by nonuniform level of the 
stray light illumination and then reconstructed the dark and 
flat-field information for individual pixels on the basis of all 
available frames. The final steps were standard dark and 
flat-field corrections. This approach resulted in a consider- 
able improvement of the photometric quality of the data. 
The implementation used our own scripts written in the 
IDL software environment. Aperture pho tometry was m ade 
by means of DAOPHOT II procedures (|Stetsonl Il987l ). as 
distributed by the IDL-astro librarjQ. 

In spite of the above careful reductions, we still observed 
linear correlations between the star flux and the sky back- 
ground level, most probably caused by a small photomet- 
ric nonlinearity of the electronic system. The correlations 
showed a trend with time which could be approximated by 
simple linear functions of time. Corrections for the correla- 
tions produced a smooth light curve of II Peg with formal 
scatter of about 0.002 - 0.004 mag. However, the light curve 
may contain slow (10 days or longer), smooth, systematic 
trends at a level of about 0.01 magnitude which cannot be 
characterized and eliminated using the available data. 

The nearby, constant, simultaneously observed stars in 
the unvigneted region of the CCD, CSC 02258-01385 and 
CSC 02258-01152, and the low amplitude S Scuti-type star 
CSC 02258-00981, discovered by MOST, were used to de- 
termine transformations between the MOST and Johnson 
V magnitude systems. The Vt and Bt magnitudes taken 

1 http:/ /idlastro. gsfc.nasa.gov/contents. html 
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Figure 2. Enlargements of eight light-curve segments with eleven flares of II Peg. Groups of data points within each figure correspond 
to individual MOST orbits. 



from the TYCHO-2 catalogue were used, after conversion 
to the standard Johnson BV system. The maximum bright- 
ness magnitude of II Peg during the first half of the MOST 
observations was estimated at V ma x ~ 7.45 ± 0.02 (random) 
with the additional uncertainty of the system transforma- 
tion of ±0.03 mag. We estimate that the combined uncer- 
tainty of the maximum F-magnitude of II Peg during the 
MOST observations does not ex ceed ±0.06. We note tha t 
the unspotted model prediction o f lBerdvueina et al l (|l998al ) 
was appreciably brighter, V u = 6.9. 

We present the light curve of II Peg in Figure [TJ where 
the V magnitudes are as determined above while the or- 
bital phases were calculated by means of the ephemeris de- 
termined by iBerdvugina et al.l (|l998al ). as quoted in Sec- 
tion Q] The accumulated uncertainty of the orbital period 
over E — 765 — 769 epochs between the original determi- 
nation and the MOST observations results in a very small 
uncertainty of the phase, ±0.002, which can be neglected in 
the present context. 

We note that during the MOST observations, the up- 
per envelope of the light curve corresponding to the Vmax 
level slowly decreased from 7.45 to 7.46, while the ampli- 
tude of light changes, AV, decreased from 0.145 to 0.12 
magnitude. It is interesting to note that the light curve ob- 
tained by MOST is similar in its shape, m aximum level an d 
amplitude to the light curve obtained by Kaluzny (1984): 
\/ =7.46, AV=0.12, as presen ted bv lBvrne etahl l| 19891 ) 
and iMohin fc Raveendranl l| 19931 ). 



(3) iMathioudakis et all (|l992h detected ten flares during 
57.4 h of optical monitoring in the Johnson U and B fil- 
ters, firidiiigJJierateof one flare per 5.9 h; 
(4a) iDovle et all (|l993l ) observed two flares in their ultravi- 
olet spectra, with one lasting about 3 hours; 
(4b) The same work reported three optical flares, lasting 
10.52, 101.00, and 9.08 min, with amplitudes 0.066 (B- 
band), 0.371 (!7-band) and 0.207 ([/-band) magnitude, re- 
spe ctively; 

(5) IMohin fc Raveendranl (|l993l ) found one flare in their op- 
tical spectra; they summarized the results obtained by other 
authors and concluded that II Peg shows a tendency to flare 
ma inly when close to i ts minimum light; 

(6) iHenrv etU] i|l996T ) estimated a flar ing rate of one flare 
per 4 .45 h, what agrees well wi th the IMat hioudakis et al.l 
(|l992T) result. They noted that iBvrne et al.l (|l994T l moni- 
tored II Peg i n 199 2 in a U filter and found no optical flares. 
IHenrv et al ] l|l996h concluded that II Peg appears to exhibit 
lon g-term changes in t he leve l of optical flare activity; 

(7) IBerdvugina et al.l ljl999ah observed two flares in optical 
spectra, with rise times of a few hours, and very long decline 
times of 1.5 and 3 d. From H n emission line profiles, they 
estimated that the flares had taken place above the visible 
pole, probably in connection with a large, single active re- 
gion; 

(8) iFrasca et al.l l)200ct l found a strong flare in spectra ob- 
tained close to light minimum, with duration time of at least 
2 d. 



3 FLARES 

3.1 Previous observations 

The astronomical literature contains a few previous reports 
of several very different flares observed for II Peg, including 
cas es of non-detect i on even for long monitoring intervals: 
(1) IBopp fc Noahl (|l980al ) observed sudden H Q enhance- 
me nts which slowly d ecayed on time scales of days; 
f2) lDovleet aD (|l99lh simultaneously detected a flare in X- 
rays and the Johnson U filter - the latter had a duration of 
more than 36 min; 



3.2 MOST results 

Three types of flares (Fig. [2} were detected during the 
MOST observations: 

- Four "short" flares (nos. 1, 2, 9, 10) lasting about one or 
two MOST orbits (2 - 4 h), and with an amplitude of about 
0.01 mag; 

- Six "long" flares, very similar in shape, rise and decay 
time, and with amplitudes of about 0.04 mag; these flares 
were used to form a "mean flare" (see below); 

- One particularly long-lasting flare, with duration of about 
one full day. 
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Figure 3. The mean flare made from five similar, long-lasting 
flares (numbers: 3, 4, 7, 8, 11), expressed in flux units. Time zero 
and maximum (1.037 continuum flux units) correspond to the 
estimated moment and mean maximum amplitude of the flare. 
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Figure 4. The long-lasting flare no. 5 together with the subse- 
quent flare no. 6 . 



Because of the high-background data gaps at 103 min inter- 
vals and the typical MOST orbit coverage of 25 min, none 
of the ten flares of the first two types was observed from 
start to end. In particular, the four short flares could not be 
analyzed sufficiently thoroughly. 

We attempted to construct a "mean flare" (in flux units) 
from the six, apparently more commonly occurring, long- 
duration flares, as shown in Figure [3] Because flare no. 6 is 
affected by the preceding unusual flare no. 5, we used the 
remaining five flares (i.e. nos. 3, 4, 7, 8, 11) for the con- 
struction. First, we expressed their intensities in contiunuum 
flux units (as denned by the underlying slow light variations 
caused by rotation of the spotted star, removed by dividing 
the data by low-order polynomials fitted to the quiescent 
parts of the light curve) and then we matched the individ- 
ual flare start times manually to an uncertainty of about 
±2 min. Then all flares were simply plotted together with- 
out any further scaling; the partially observed flares nos. 3 
and 8 contributed only the decaying parts to such a mean 
flare. 

The estimated rise time from the flat continuum to the 



0.86 



phase 0.9-2.1 
phase 3.9-5.1 




0.1 















0.2 


).3 


0.4 0.5 


O.f 


i 0.7 


0.8 



phase 



0.9 1 



Figure 5. Individual flares are indicated by vertical marks along 
the horizontal (orbital phase) axis. Two sections of the light curve 
of II Peg, from the beginning (orbital phases 0.9 - 2.1) and the 
end (phases 3.9 - 5.1) of the MOST monitoring are shown for 
comparison. The light curve evolved mainly due to a drift in lon- 
gitude of the spots at a rate faster than the orbital motion. 



maximum of the mean flare was found to be 25 ± 4 min. The 
decline time, from the maximum back to the flat continuum, 
varied in the range of 5 to 10 h. The half-maximum duration 
time of the mean flare was one hour (To. 5 « 59 min, as de- 
fined in iKunkell l|l973h ) spanning the range between 50 min 
for flare no. 8 and 74 min for flare no.7. The duration time 
was shorter for flare no. 11, but it could not be uniquely de- 
termined from the available data. All these flares considered 
here most probably did not occur on the secondary com- 
ponent of II Peg, the M-dwarf, for which values o f T0.5 of 
the o rder of hundreds of seconds would be expected ()Kunkell 
1973h : the location was most likely the primary component 
or somewhere in the space between the stars. 

In general, the first two types of flares observed by 
MOST , were similar to those observed before bv lDovle et al.l 
|l993h (items (4a) and (4b) in the previous section). 
The long-lasting flare (no. 5), which started at HJD — 
2, 454, 738.6, close to light minimum, may be an analogue of 
flares ob served to date as enha n cements of optical s pectra l 



Hares observed to date as enhancements ol optical spectral 
lines by Bopp fc Noahl l|l980al ). iBerdvueina et al.l (|l999af ) 
and iFrasca et all (|2008r i. It differs markedly in its shape, 
rise time (6 h) and decay time (at least 18 h, possibly 24 h) 
from the remaining ten flares observed by MOST. It is shown 
among the other flares in Figure[2]and magnified in Figure[4] 
The rate of eleven flares in the time span of 30.877 d 
gives a flaring rate of about one flare per 2.8 d. However, due 
to the breaks in the MOST observations, we cannot neglect 
the possibility of overlooking very short flares lasting only a 
few minutes ; these would be fla res similar to the two shortest 
observed bvlDovle et al.l (1 19931) and all of those observed by 
iMathioudakis et all (|l992T ). 



3.3 The phase distribution of flares 

The phase distribution of flares in relation to the spot- 
modulated light curve can be inspected in detail in Figure [5] 
The flares appear not to be uniformly distributed in orbital 
phase: as many as five flares appeared within the light mini- 
mum, in the orbital phase interval 0.75 - 0.85. This supports 
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Figure 6. Synthetic light cur ves for II Peg computed u sing the 
stellar parameters obtained bvfecrdvugin a et al. (1998a) for four 
different values of the orbital inclination. 



the conclusion of lMohin fc Raveendranl (| 1993h that flares in 
II Peg are concentrated close to light minimum when the 
most heavily spotted side of the visible star is directed to- 
ward the observer. However, a Kolmogorov-Smirnov test for 
the deviation of the phase distribution from uniformity gave 
the probability that the distributions appear to be identical 
of 0.28. While this is a small number, it is not small enough 
to prove this assertion, which still requires confirmation. 



4 A SEARCH FOR ECLIPSES 

High-precision photometry from space carries a potential for 
detection of eclipses caused by transits of the undetected 
secondary companion over the visible star. We estimated 
the expected depths and durations of primary eclipse in the 
Johnson V-band for several values of inclina tion using the 
Wilson-Devinney light curve synthetic c ode ([Wilson! Il996h 
for the physical parameters obtained by iBerdvugina et all 
(|l998al ); this is shown in Fig. |6] and Fig. [7] 

A careful inspection of the MOST data revealed no in- 
dication of any eclipses, as can be seen in Figure To an- 
alyze the conjunction segments of the light curve, trends 
introduced by spots were fitted within phase ranges 0.93 - 
0.97 and 1.03 - 1.07 and then normalized light curves were 
analyzed in great detail. The data do not reveal any system- 
atic, localized deviations which would have depths similar to 
those predicted in Fig. [5] to less than 0.1 per cent. 



5 LIGHT CURVE ANALYSIS 

In this investigation we modeled the MOST light curve in 
terms of dark spots which are internally invariable in time. 
We used the pr ogram StarSpotZj which was su ccessfully ap- 
plied to e Eri (|Croll et al.ll2006l ) and k 1 Cet (jWalker et all 
120071 ). where differential rotation of the stars was found. 
The reader is directed to these papers for details of the 
model. The p rogram is based on the program SpotModel 
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The assumption of the internal invariability of spots 
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Figure 7. Segments of the MOST data close to predicted times of 
primary eclipse. The segments for four consecutive conjunctions 
(which show spot evolution in time) have been shifted down for 
clarity by the indicated magnitude amounts; the first conjunction 
is on top. The start (90s, 77s) and end (77e, 90e) phases of the 
eclipse, estimated for i = 90° and i = 77° are represented by four 
vertical lines. Note that, as discussed in Sec. [2] the current spec- 
troscopic ephemeris is very accurate and predicts the conjunction 
phase to ±0.002. 



is most likely not fulfilled in the case of I I Peg. As Doppler 
images obtain ed during 1994-2002 reveal (|Berdvugina et all 
Il998bl . fl999bh . the spots may constantly change their prop- 
erties over time, in time scales of several rotation periods. 
The shortest time scales for appreciable spot changes could 
be as short as 2 months, which is comparable with the length 
of the MOST run of 31 days. The global, progressive light 
curve shape changes (Fig. [TJ can be easily explained by dif- 
ferential rotation of the stellar surface with small, random 
changes of the spots well averaged over the time of observa- 
tions. 

A more detailed investigation of differential rotation should 
be supported by simultaneous high-resolution spectroscopic 
observations. Without them, as in the current investigation, 
we are unable to assert whether spots really remained suffi- 
ciently constant during the MOST observations. Addition- 
ally, whatever method of spot shape restoration is used (in- 
cluding the maximum entropy method), the results will al- 
ways be subject to limitations imposed by the instrumental 
effects mentioned in Section [2] 



5.1 Light curve modelling 

After first trial runs it turned out that at least two cold spots 
on the hemisphere directed to the observer are necessary to 
give a reasonable explanation of the light variations. 
Compatible with the non-detection of eclipses, we assumed: 
i = 60°, Ri = 3.4 R R and v sint = 22.6 km/s, as de- 
termined by IBerdvugina et all (|l998a) (see also fig. 6 of 
their paper). Also, as fixed parameters for the light curve 
models, we assumed the linear limb- darkening coefficient 
u = .817, adopted using the tables of iDfaz-Cordoves et alj 
(1995). We also fixed the photospheric and spot tempera- 
tures at T phot = 4, 600 K and T a p t = 3 , 600 K , as based 
on the results of IBerdvugina et al. I (|l998bl.[l999bl ). This as- 
sumption was required to fix the spot-to-photosphere flux 
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Table 1. Results of the light curve models for a rigidly and a differentially rotating star for i = 60°. The ranges in resulting 
parameters for the two extreme values of the normalization parameter F u , 1.19 and 1.33 (see the text) are given in brackets as 
estimates of parameter uncertainties, 
a - assumed as constant during modelling, 

b - determined using the constraints Ri = 3.4 Kq and i;sini = 22.6 km/s, 
c - calculated using Eq.(l). 
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- 28.645) 


28.166 (28.045 - 28.220) 


4>2 [°] 


18.1 (13.6 - 


21.3) 


75.7 (55.0 - 79.7) 


21.2 (17.7- 


24.4) 


71.7 (49.1 - 76.6) 


T2 [°] 


15.0 (14.6 - 


15.4) 


27.4 (15.5 - 36.3) 


16.5 (16.0 - 


17.2) 


26.4 (15.6 - 34.9) 


PS [d] 














t 3 [hjd] 














03 [°] 


-90 a 




-90 a 


-90 a 




-90 a 


ra [°] 


89 a 




89 a 


89 a 




89 a 


Y 2 , 

A- red. we i gh 


13.79 




11.75 


14.72 




12.28 



ratio / for the bandpass of the MOST observations. / = 
0.077 ± 0.015 has b een evaluate d by means of the SPEC- 



TRUM programme ( Gray||200ll ) and Kurucz's atmosphere 
models (|Kurucdll993l ) using the MOST filter bandpass. The 



techni q ue: According to iNeff et alj ill995l). lO'Neal fc Neflj 
l|l997t ). [Marino et all (l999h and lO'Neal et all (|l998l ) spots 



same value of / was assumed for all spots. The remaining 
parameters of the model were, for each spot: (1) the initial 
moment t (in hjd = HJD — 2, 454, 700), when the spot is 
exactly facing the observer, (2) the rotation period of the 
spot p in days, (3) the latitude <f), (4) the diameter r in de- 
grees and (5) the value of unspotted flux F u — 1.26. 
The latter value was calculatted assuming the value of 
unspotted magni tude equal to V u = 7.20, as determined by 
IChueainovl l|l976t ) at the time when he observed a flat max- 
imum, and a lso adopted by Mohin fc Raveendranl l)l993h . 
We note that lBerdvugina et al. ( 1998bl ) suggested V u = 6.9, 
but for such a high brightness it is impossible to obtain a 
physically plausible fit to the light curve: We would have to 
postulate that we observed II Peg with almost the whole sur- 
face covered by black spots. As we described in Section[2] we 
observed V ma x = 7.45±0.06, based on the MOST instrumen- 
tal system, after transformation to the V-band using nearby 
stars. If the unspotted magnitude for II Peg is V u = 7.20, the 
unspotted flux at the time of the MOST observations would 
be F u — 1.26 ± 0.07 (using normalization of F u = 1 for 
V m ax = 7.45). This leads however to the difficulty of large 
radii of both spots with some overlap, which is not admit- 
ted by the model. We solved this problem assuming a third, 
very large (r$ — 89°) circular spot, covering practically the 
whole hemisphere directed away from the observer. This 
spot remained constant and - because of the low inclina- 
tion - only partially visible. It represented the non-variable 
part of the spotted photosphere. This assumption is strongly 
supported by results obtained from the Doppler imaging 



are always visible and they cover between 35 to 64 per cent 
of the hemisphere projected toward the observer. As men- 
tioned in Section [2] the amplitude of light changes observed 
by MOST was AV = 0.145 - 0.12 magnitude. This was 
close to the smallest value noticed to date which, accord- 
ing to iMohin fc Raveendranl l|l993h . means that during the 
MOST observations spots covered a large fraction of the 
stellar surface. 

The light curve used in the model utilized mean points 
formed for individual MOST orbits, after removal of all stel- 
lar flares. Typically 40 - 70 points contributed to one MOST- 
orbit average point. The formal normalized flux errors (a) 
per point is about 0.0011 (median) and the full range of 
0.0007 - 0.0022. Please note, that as discussed in Section [2] 
the light curve may contain a roughly 10 day-long smooth 
trend, a few times larger in amplitude than the formal nor- 
malized flux errors given above. 

The spot model using the StarSpotz program assumes 
spherical stars. II Peg is a binary system where some light 
modulation is expected from ellipsoidal and reflection ef- 
fects, so-called proximity effects. They are small in an abso- 
lute sense (less than 0.02 mag, see Figure [5} but cannot be 
neglected as, for the assumed i — 60°, they reach almost 9 
per cent (0.012 mag) of the observed spot-modulation ampli- 
tude. Because of the uncertainty with the inclination (±10°), 
all physical parameters of II Peg are not fully known; this 
affects predictions of the proximity effects. To assess the im- 
pact of the proximity effects on the final results, particularly 
on the differential rotation of the visible star, we used two 
light curves, corrected and uncorrected for the proximity ef- 
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fects. Each of the two light curves was analyzed assuming a 
rigidly and a differentially rotating stellar surface. 

The individual spot rotation periods pi, 2 were assumed 
to depend on the stellar latitude (j>i,2, the rotational pe- 
riod on the equator P eq and differential rotation coefficient 
k through: 

Pi{4>i) — Peq I {I - fcsin 2 4>i), (f) 
where i=f,2. 

The search procedure for k consisted of two steps: first, 
for the assumed R — 3.4 R©, the proper value of P eq = 
6.5940±0.0005 d was found, providing the observed v sin i = 
22.6 km/s; then - for the value of unspotted flux level 
F u = f .26 - the differential rotation coefficients k, return- 
ing the smallest value of the reduced and weighted \ 2 was 
derived. The formula given by Eq.(l) corresponds to the 
solar-type differential rotation law. Due to the low quality 
of our fits, most probably dominated by the inadequacies of 
the spot model (see Section f5.2p . we did not consider other 
possible types of differential rotation. 

5.2 Results of the light curve modelling 

The results of modelling are presented in Table Q] One can 
see in Figure [8] that the fit obtained for the case of differen- 
tial rotation describes the light curve better than the solid- 
body rotation model. This applies to the general light curve 
evolution, and in particular to the progressive amplitude de- 
crease, as discussed in Section [2] We also note that only for 
the differential-rotation models does the l arger spot face the 
secondary star , similarly as obtained by iBerdvugina et all 
l|l998bl . fl999bh . 

To estimate the systematic errors of our models result- 
ing from the large uncertainty of V m ax of ±0.06 mag, we 
repeated our solutions for two values of the unspotted flux 
level, F u = 1.19 and F u — 1.33. This choice affects the solu- 
tions strongly and the resulting spread in parameters can be 
taken as an indication of the uncertainty in our solutions. 
Note that the order in the range limits given in Table Q] is 
sometimes inverted but the first value always corresponds 
to the smaller value of F u . 

In general, the residuals - typically at a level of 0.004 
of the mean flux - are much larger than formal errors of in- 
dividual data points (typically 0.001). This has driven val- 
ues of the formally derived, reduced, weig hted x 2 (Table fTJ) 
to values well above unity indicating systematic trends in 
residuals, most probably reflecting the difference between 
the true and the circular shape of spots which was assumed 
in the model. 

Because of the dominance of the systematic deviations over 
the random noise in the values of x 2 , this parameter has 
only an indicative utility. Nevertheless, for each pair of so- 
lutions, with included (columns 3, 4) or excluded (columns 
1, 2) proximity effects, the differential-rotation solution ap- 
pears to be always better than the solid-rotation one. Taking 
these considerations into account, we select the solution in 
the last column of Table [1] as the final one, and we plot it in 
Figure [8] 

5.3 Comparison with other results 

iHenrv et al.1 j 19951 ) determined the differential rotation pa- 



rameter, k = 0.005 ± 0.001, for II Peg using several multi- 
epoch light curves. Our result of k = 0.02451q qJ^o i s m 
better accordance with the linear rel ation between param- 
eters of RS CVn-type stars (Eq. 9 in IHenrv et al.l (|l995h ^ : 
log k = -2.12(12) + 0.76(6) x log P rot - 0.57(16) x F, where 
F = Rstar / RrocHb- Using the parameters listed in Table [TJ 
we have P rot ~ 6.7 d, F = 3.4/7.1 « 0.48, leading to a pre- 
diction of k — 0.017. However , when we take int o account 
the scatter visible in fig. 28 of IHenrv et al. I l|l995h . a broad 
range of 0.002 < k < 0.066 is admitted for this value of 
Prot- Interestingly, the value of k determined in this paper 
for II Peg is similar to that estimated for the apparently 
single, but even f aster rotating giant, F K Com (k = 0.016 
for P = 2.4 d) bv lKorhonen et al.1 l|2002i ). 



6 CONCLUSIONS 

Analysis of the almost-continuous, one month-long photo- 
metric monitoring of II Pegasi by the MOST satellite per- 
mits us to formulate the following conclusions: 

(i) Eleven flares were observed, one lasting about 24 h and 
six flares moderately long, lasting typically 5 to 10 hours. 
The characteristics of the four shortest flares were difficult 
to estimate. 

(ii) The primary eclipse of the visible star by its compan- 
ion (probably M-dwarf) was not detected, which gives an 
upper limit for the orbital inclination of the system of 76° . 

(iii) From the analysis of the dark-spot modulated light 
curve, assuming i = 60°, R\ = 3.4 R0, vsini — 22.6 km/s 
|Berdvugina et aTlll998al ) and absence of internal variability 
of spots during the MOST observations, we obtained an es- 
timate of the parameter measuring the differential rotation 
of the primary component of II Peg: k = 0.0245±o;°j2o • The 
error of k reflects the major uncertainty in the unspotted 
brightness of the star so that the value of k remains prelim- 
inary; it will improve with future ameliorations in values of 
the assumed stellar parameters which enter the model. 
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Figure 8. The fit to the light curve of II Peg (corrected for proximity effects) by the model with solid-body rotation (k = 0, column 3 
of Table \Tl and with differential rotation (k = 0.0245, column 4). 
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